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The non-classical nature of quantum states, often illustrated using entanglement measures or 
quantum discord, constitutes a resource for quantum information protocols. However, the non- 
classicality of a quantum system cannot be seen as a property of the state alone, as the set of available 
measurements used to extract information on the system is typically restricted. In this work we 
study how the non-classicality of quantum measurements, quantified via their incompatibility, is 
influenced by quantum noise and how a non-Markovian environment can be useful for maintaining 
the measurement resources. 
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I. INTRODUCTION 

The non-classical nature of quantum states is consid¬ 
ered to be an essential resource for the emerging quantum 
technologies. This idea is well-established in the context 
of entanglement theory |l], and has also been formulated 
generally in the framework of quantum resource theo¬ 
ries 01 . However, characterising the non-classicality of 
a quantum system as a property of the state alone has 
limited practical significance, since the set of available 
measurements used to extract information on the system 
is almost always restricted by experimental limitations. 
Accordingly, measurement resources are expected to play 
a significant role in realistic quantum devices. 

The study of measurement resources is motivated also 
from a fundamental point of view. As shown in [4| there 
exist entangled states which nevertheless cannot be used 
to produce quantum correlations in any Bell experiment 
because of the existence of a specific hidden variable 
model for the correlations. For certain scenarios based 
on correlation experiments, the existence of such a clas¬ 
sical description is equivalent to joint measurability of 
suitable observables. In particular, it has been recently 
shown that for the CHSH-Bell inequality [5] and EPR- 
steering scenarios |6-[8j \W \, the appropriate quantum re¬ 
source can be formulated in terms of incompatibility of 
the available local observables. It is therefore important 
to aim for a better understanding of the incompatibility 
of observables. 

Since every real quantum system interacts with its en¬ 
vironment mm , a practical implementation of any 
quantum protocol has to take the possible effects of noise 
into account. The dynamical behaviour of nonclassical 
features of quantum states has been extensively studied 
in recent years. Specifically, understanding dynamical 
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phenomena such as entanglement sudden death [T3], or 
frozen discord m has been useful for determining the 
level of isolation from the environment required to im¬ 
plement protocols relying on these quantum features. 

In order to investigate the dynamical loss of measure¬ 
ment resources, we need an approach somewhat different 
from the conventional state-centred view: decoherence 
and dissipation now take place in the Heisenberg picture , 
rendering the measurements less useful for revealing non- 
classical properties. In fact, consider the general setting 
consisting of preparation of an entangled state shared by 
Alice and Bob, followed by noisy local evolution (quan¬ 
tum channel) on Alice’s side, and local measurements 
performed by Alice and Bob at the end. In this setting, 
one would usually write the associated joint probabilities 
in the Schrodinger picture, with the local noise induc¬ 
ing loss of entanglement on the state. However, we can 
equivalently use the Heisenberg picture, acting on Alice’s 
measurements. 

An essential starting point for our investigation is the 
fact that incompatibility can never be created but is of¬ 
ten destroyed by the action of a quantum channel m- 
How exactly this “measurement decoherence” takes place 
during actual open system dynamics has not been inves¬ 
tigated before and the purpose of this paper is to take 
the first step in this direction by performing a quanti¬ 
tative analysis of the dynamical evolution of incompat¬ 
ibility. We concentrate on two well known microscopic 
open system models, where the dynamics can be tai¬ 
lored via environment engineering techniques: the phase 
damping m and amplitude damping El evolutions un¬ 
der both Markovian and non-Markovian noise. We find 
that both evolutions exhibit sudden death of incompat¬ 
ibility of quantum measurements, even though no en¬ 
tanglement sudden death occurs. We further study the 
case of a highly engineered amplitude damping dynam¬ 
ics, given by a photonic band gap environment, which 
has been found to efficiently protect entanglement in the 
long time limit El- We find that even for such a highly 
engineered scenario the incompatibility of relevant mea¬ 
surements cannot be maintained. This demonstrates that 
incompatibility is more fragile than entanglement, which 
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is in line with the recently observed fact that an entangle¬ 
ment breaking channel is incompatibility breaking, but 
not vice versa HE. m- 

The structure of the paper is as follows. In Sec. [11 
we review general features of quantum incompatibility, 
as well as the specific quantification introduced in m- 
In Sec. [HI] we motivate the following sections by a simple 
classical noise model. In Sec. [IV]we proceed to describe 
incompatibility under quantum dynamical noise, and in 
Sec. [V]we introduce the microscopic open system models 
under consideration. The numerical results are presented 
and discussed in Sec. |vg and we conclude the paper in 
Sec. EII1 


II. QUANTUM INCOMPATIBILITY 


probability distribution p = (jpi). The POVM E A,P de¬ 
scribing the resulting deformed measurement is then 

E?’ p = (1 - A )Ei + X Pi t, (1) 

and for a given pair (E, F ) of measurements, we look for 
the minimal value of A for which the noise-deformed pair 
(E A,P , F x,p ) becomes compatible: 

Z p (E, E) := inf{A > 0 | (E A,P , F x,p ) compatible }. (2) 

It is easy to see that this function fulfils the properties 
(i) and (ii) of an incompatibility monotone, and we also 
know that A < 1/2 [24 . There are also other choices for 
incompatibility monotones, of which we have included a 
brief discussion in Sec. VII. 


A. Definition of incompatibility and its 
quantification 

We first recall that a general quantum measurement is 
described by a POVM, which (assuming a finite set of 
measurement outcomes) is a map E : i E{ assigning a 
positive operator Ei to each outcome i and satisfying the 
normalisation condition JV E^ = 1. Two measurements 
E and F are compatible or jointly measurable if there ex¬ 
ists a measurement G = ( Gij ) such that Ei = JV Gij 
and Fj = JV Gij . Intuitively, this means that the mea¬ 
surements can be realised together with a single device. 
The definition extends naturally to collections contain¬ 
ing more than two measurements. For the purposes of 
the present investigation, however, studying the evolu¬ 
tion pairs is sufficient. If the measurements are not com¬ 
patible, they are said to be incompatible. In this case 
they can potentially be used for producing nonclassical 
features as discussed in the introduction. 

In order to quantify incompatibility, some of the au¬ 
thors of this paper recently introduced m the concept 
of a incompatibility monotone as a function X on pairs of 
observables, required to satisfy the following conditions: 

(i) Z(E, F) = 0 if and only if E and F are compatible. 

(ii) X(A (E), A(E)) < Z(E, E) if A is any channel. 

The channel A in (ii) can be either a quantum channel 
applied before the measurement, or a classical channel 
applied on the measurement outcomes. This definition 
of an incompatibility monotone is analogous to entangle¬ 
ment monotones pQ. 

In order to define an explicit quantification satisfy¬ 
ing the above conditions, one can use m the general 
resource-theoretic notion of noise-robustness [21 f20H24] . 
Basically the idea is to use classical selection noise as 
a reference. In our case, this can be operationally im¬ 
plemented as follows: we think of implementing a mea¬ 
surement E using a noisy device which sometimes (with 
probability A) ignores the actual outcome, and instead 
outputs randomly an outcome i according to some fixed 


B. Incompatibility in a qubit system 

In this paper we focus exclusively on single qubit sys¬ 
tems for which the optimisation <© can be reduced to 
solving a single polynomial equation (see discussion in 
m)- In order to describe this, we first need to recall the 
representation of qubit measurements. 

Binary POVMs have only two outcomes, 0 and 1. Due 
to the normalisation, we have Eo = 1 — Ei, so that the 
operator 0 < E\ < 1 in fact specifies the POVM com¬ 
pletely. In a qubit system we can then write E\ in terms 
of its Bloch four-vector x = (xq, aq, aq? x s) = (a?ce x ) : 

Ei = ^(x 0 t + x • a), (3) 

where a = (01,02, 03) and cr/s are the usual Pauli ma¬ 
trices. The Bloch vector for the other POVM element Eg 
is given by x 1 - := (2 — a?o, —x). The condition 0 < E\ < 1 
can be compactly expressed by r, r 1 G E+, where 

E + = {x | (x\x) > 0, xq > 0} (4) 

and 
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(x\y) := x 0 y 0 - Y x iVi ( 5 ) 

is the Minkowski scalar product. To summarise, a binary 
measurement in a qubit system is specified by a single 
four-vector x satisfying x^x 1 - G T+. 

In general, deciding whether a pair of measurements is 
compatible requires solving a convex optimisation prob¬ 
lem. In a qubit system, the problem simplifies consider¬ 
ably [25] (see also pa mi) : two measurements, specified 
by Bloch vectors x and y, are compatible if and only if 

c (x,y) > 0, 

where 

c (x,y) := [{x\x){x ± \x- L )(y\y){y ± \y ± )] 1/2 - {x\x ± ){y\y^) 
+ (x\y J -){x- L \y) + (x\y){x ± \y ± ). 
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Given any two measurements with Bloch vectors x and 
2 /, we can now write down the formula for the noise- 
robustness-based quantification § of their incompati¬ 
bility. Since we only have two outcomes, the classical 
noise distribution p is written conveniently in terms of 
the bias parameter — 1 < b < 1: 

Pi = ^a + b), p 0 = ^(l-b). (6) 

Hence, the noise-deformation 0 transforms a given 
Bloch vector x as 

x N x ,b(x) := ((1 - X)x 0 + A(1 + 6), (1 - A)x), 

so the incompatibility quantification for a pair (x,y) is 
given by 

Ib(x,y) ■= inf {A > 0 | C(N x ,b(x), N x , b {y)) > 0} • (7) 

Since Xb(x,y) < 1/2, and given that x and y are incom¬ 
patible, we know that A C(N\^(x), N\^(y)) changes 
sign at exactly one point on the interval [0,1/2], and this 
point is Zfc(£, y). This value can be easily determined nu¬ 
merically for any given pair x, y, as a root of a polynomial 
equation. 

We close this section by recalling an alternative oper¬ 
ational interpretation of the number Z&(x, y) for the case 
b = 0. In fact, it coincides with the maximal violation 
of the CHSH-Bell inequality achievable with Alice’s mea¬ 
surements fixed to be x and y Ha¬ 
ul. INCOMPATIBILITY UNDER CLASSICAL 
NOISE 

In order to illustrate the relationship between entan¬ 
glement and incompatibility in noisy systems, it is help¬ 
ful to begin with a simple case where the noise is purely 
classical, obtained by mixing with a completely depolar¬ 
ising noise. This is a variant of the selection noise used 
already in the preceding section to define the incompat¬ 
ibility quantification, and appears in the construction of 
the well-known isotropic states pQ 

Q s = (l-s)\* 0 ){* 0 \ + sl/d 2 , 0 < s < 1, (8) 

where T 0 is a maximally entangled state of a bipartite 
system. The parameter t describes classical random se¬ 
lection noise added to the maximally entangled state, 
and one can show that the state loses its entanglement 
at some critical value s = ^entangled? remaining separa¬ 
ble for s > s en tangied- However, the state loses its non- 
classicality already before , in the sense that after some 
S = Squantum < ^entangled a hidden state model [1 HD] 
exists for any bipartite correlations of the form 

W s (i,j\a, b) = tr [g s A“ 0 B b ], (9) 

where A 1 , ..., A n and B 1 ,..., B 171 are local measure¬ 
ments accessible to the two parties Alice and Bob, re¬ 
spectively. Hidden state models are specific hidden vari¬ 
able models for which the correlations P s (i,j|a, b) do not 


violate any Bell inequality. Moreover, the nonexistence 
of such a model is equivalent to the possibility of EPR- 
steering HD], which has practical applications to e.g. one¬ 
sided secure Quantum Key Distribution (QKD) and sub¬ 
channel discrimination [28] [29]. 

For a fixed state, the quantum resource for steering 
can be entirely associated with measurements. In order 
to see this, we first write 

0»=(A:®l)(|¥ o )<¥o|), (10) 

where the effect of the classical noise is described by a 
local quantum channel A s = (1 — s)Id + st/d and A* 
denotes the Schrdinger picture. Now, since the noise and 
measurements are both local, it makes more sense to use 
the Heisenberg picture: 

P s(i,j\a,b) = ('l>o\K(A?)®B b j * 0 )=tr[A s (A?) T B b j }. 

( 11 ) 

Here, the state Tq is maximally entangled, the transpose 
is taken in the associated basis, and we observe that the 
effect of noise is entirely captured by its action on Alice’s 
measurements. One should notice the following impor¬ 
tant observation [6j El E] : the nonexistence of the hidden 
state model for the correlations P(i, j\a, b) is equivalent to 
incompatibility of Alice’s noisy local observables 

A s (A a ), a = 1 ,..., n. ( 12 ) 

In particular, this means that the observables A s (A a ) lose 
their incompatibility at the amount of noise s = s quantU m 
where the hidden state model comes into existence. For 
s < ^quantum 5 the measurements have some degree of in¬ 
compatibility, which can be quantified as described in 
the preceding section. In the next section, we proceed to 
investigate this loss of incompatibility under more com¬ 
plicated quantum noise channels arising from dynamical 
interaction with a quantum environment. 

IV. INCOMPATIBILITY UNDER QUANTUM 

DYNAMICAL NOISE 

Evolution of an open quantum system can be described 
by a family (At) of quantum channels indexed by the 
time parameter t > 0. This map is typically obtained 
by assuming unitary dynamics t i— U t on a total system 
consisting also of an environment, initialised in a given 
state <r env . In fact, given that the system is initially (time 
zero) prepared in a state g, and measured at time t via a 
POVM (Ej), the probability for the outcome j is given 
by 

ti[gA t (Ej)] ■= tr[Ut(g® ^env Wt(Ej 0 lenv)] • (13) 

This determines a family of completely positive unital 
maps A t acting on the observable algebra of the system, 
describing the evolution in the Heisenberg picture. That 
is, (A t{Ej)) is the POVM describing the measurement 
(Ej) performed at time £, from the point of view of the 
initial state g. 
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A. Dynamics of incompatibility 

The properties of the measurements may crucially 
change depending on the point in time at which the 
measurements are performed. In particular, incompat¬ 
ibility, as quantified by any incompatibility monotone X, 
becomes a function of time: 

t^I(A t (E),A t (F)). (14) 


so that 

xo(t) = x 0 + (a(t) - l)x 3 
x\(t) = Re[c(t)]xi — lm[c(t)]x 2 
x 2 (t) = lm[c(t)]xi + R e[c(t)]x 2 
x 3 (t) = a(t)x 3 . 

Given two initial vectors x and y , the incompatibility 
evolution 


As a motivating starting point for our investigation, 
a few general remarks can be made concerning the case 
where t A t is a continuous semigroup , i.e., 


At+t'— At(A t '(-)). (15) 

In this case it follows directly from the general properties 
(i) and (ii) of an incompatibility monotone that the func¬ 
tion ( pA[ ) is nonincreasing. Furthermore (assuming finite 
system dimension), in a generic (ergodic) case there ex¬ 
ists a unique stationary state g s , such that 

lim A t (-) = tr[^ s (-)]l, (16) 

t —yoo 


t I b (x(t),y(t)) 


( 20 ) 


can then be computed as explained in Subsec. |IIB| 
Throughout the calculations, if not mentioned other¬ 
wise, the initial measurement pair is chosen as the spe¬ 
cific maximally incompatible pair with xo = yo = 1, 
x\ = 2/2 = 1 , x 2 = yi = 0 and x 3 = y 3 = 0 (denoted 
Pi). We also study the case xq = yo = 1, x\ = y 3 = 1, 
x 2 = U2 = 0 and x 3 = yi = 0 (denoted V 2 )- Further¬ 
more, we have set the bias parameter 6 = 0. More general 
measurements could be considered in a similar fashion. 
However, these two cases already exhibit several interest¬ 
ing phenomena, as we will see below. 


which implies that at the limit, every measurement is 
trivial. Hence one expects that the incompatibility func¬ 
tion (14) will eventually decrease to zero under the as¬ 
sumption (15). 

The purpose of the present paper is to study the be¬ 
haviour of the function (14) for certain concrete choices 
of the dynamical map A t on two-level systems, using 
the incompatibility monotone 0. More general scenar¬ 
ios could easily be handled with more numerical effort. 
However, the two-level system already exhibits the rele¬ 
vant phenomena to the extent sufficient for the illustra¬ 
tive purposes of the present paper. We represent a given 
dynamical map t i— A t in the Heisenberg picture, and 
then compute the associated evolution t i—A t (E) for 
any given measurement T, in terms of its Bloch vector 
x(t), where x(0) specifies the initial operator E\. 

The dynamical maps considered in this paper act in 
the Schrodinger picture on a matrix T = [T^] as 


/T 0 o T 01 \ f a(t)T 0 o c(t)T 0 1 \ / 17 x 

V^io T X1 ) ^ \c(tyT 10 Tu + (1 — a(t))TooJ ’ [U) 


where a(t) and c(t) are explicit functions of time t. The 
Heisenberg representation of this evolution is determined 
by the equation 


tr[A*(T)5]=tr[TA t (5)], (18) 

where T and S are arbitrary two-by-two matrices. Fixing 
S = ^(xol+x-cr) corresponding to an initial Bloch vector 
x = (xo,x), we then find the Bloch vector x(t) of A t (S) 
via 


xo(t) = tr[A*(l) 5], Xi{t) = tr[A*(<Tj) 5], (19) 


B. Comparison with entanglement dynamics 


It was shown in section |III| that incompatibility and 
entanglement can be compared in a bipartite scenario 
involving local measurements and a shared state. More 
specifically, we look at the scenario from the point of view 
of one party, say, Alice, and compare the incompatibility 
of her measurements with the entanglement of the shared 
state. Bob’s measurements do not play a role in this set¬ 
ting. We remark that such asymmetry appears naturally 
in one-sided protocols such as steering m mentioned in 
the Introduction. 

Accordingly, we take an ancillary qubit (in addition 
to our open qubit system described above), and set the 
initial state of the combined system to be the maximally 
entangled state T 0 = ^(|00) + |11)). We stress that this 
specific choice is not crucial, since the connection between 
incompatibility and entanglement holds for any bipartite 
state [9b Analogous to eq. (10), we then consider the 
evolved state 


Qt = (A t *^Id)(|T 0 )(T 0 |), (21) 

where A^ is now the local noise channel of the form 
According to the discussion in Sec. m loss of incom¬ 
patibility of any set of Heisenberg-evolved local observ¬ 
ables can now be compared with the loss of entangle¬ 
ment of the state g t . In particular, we know from the 
general argument that incompatibility is lost before en¬ 
tanglement. In order to compare finer details of the re¬ 
spective dynamics, we need to choose some quantification 
also for entanglement of g t , and in our investigation we 
will use concurrence [30]; this has the advantage of con¬ 
necting our work with the extensive existing literature on 
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entanglement dynamics. Moreover, we would like to note 
that for the models considered here, the entanglement 
dynamics happens to be equivalent to the dynamics of 
the information flux m the non-monotonicity of which 
indicates the non-Markovian character of the dynamics. 

A straightforward calculation shows that the matrix 
of the evolved state g t has nonzero elements only in the 
diagonal and antidiagonal entries. Thus, the concurrence 
S(gt) has the simple form [32] 

£(Qt) = 2m.ax{0,K 1 (g t ),K 2 (et)}, ( 22 ) 

with 


Ki( Qt ) = | (00| Qt 111) I - y<10|ft|10> (01| ~Qt |01) 

K 2 {Qt) = I (Oil Qt |10) | - V(00| & |00) (111 Qt 111), 

leading to 


£(<h) = m |. (23) 


V. DESCRIPTION OF THE DYNAMICS 

We now describe the two different choices for the dy¬ 
namical map in detail, including examples of physically 
relevant scenarios. 


A. Phase damping dynamics 


Let us first consider a dephasing model, involving only 
decoherence without dissipation. This type of evolution 
can be experimentally implemented, for instance, in an 
optical setup [34] . where the environment can be easily 
tuned so as to produce a variety of dynamics both in the 
Markovian and non-Markovian regimes. In fact, the time 
evolution of a photon traveling in a quartz plate may be 
described by a unitary operator U ( t ) defined by 

E7(i)|A)<8> |w) = e inxut \\)® |w), (24) 


where n\ represents the refraction index for light with 
polarization A = H,V. The presence of the quartz 
plate thus leads to the decoherence of the superposi¬ 
tions of polarization states, due to the nonzero difference 
An = ny — nn in the refraction indices of horizontally 
and vertically polarized photons. The corresponding dy¬ 
namical map AJ in the Schrodinger picture takes the form 
(17) with 


a(t) = 1 , c(t) = n(t ), (25) 

and the decoherence function n(t) is given by the Fourier 
transform of the frequency distribution of the photon, 

K (t) = jdu\f(u)\ 2 e iuAnt . (26) 


The frequency distribution is modified in the experi¬ 
ment via variation of the tilting angle of an FP cavity in¬ 
serted in the path of the photon [34]. All such frequency 
distributions are very well approximated by a sum of two 
Gaussians centred at frequencies ujk with amplitudes 
and equal widths <r. This yields 

K(t) = 4_e-5 ff 2 t 2 (e-“ lt + Ae~ iu)2t ) (27) 

where A\ = — 7 ^-, A 2 = yy^ and A 00 = co >2 — The 
tilting angle of the cavity is relatively small and thus the 
distance Auj between the peaks remains approximately 
constant. 


B. Amplitude damping dynamics 


Let us now consider a case where the interaction be¬ 
tween the quantum system and its environment also ex¬ 
hibits energy exchange, that is, the interaction is dis¬ 
sipative. We look at a microscopic Hamiltonian model 
describing a two-state system interacting with a bosonic 
quantum reservoir at zero temperature given by mi 

H = tu 0 <j z + yy fe a).afe + ^(g k a k a + + gla\o-). (28) 

k k 


As usual, a± are standard raising and lowering operators 
respectively. The dynamics of a single amplitude damped 
qubit is captured by the time-local master equation: 


dQt 

dt 


= i (t) 


a -Qt&+ — Qt) , 


(29) 


where a± are the spin lowering and rising operators. The 
resulting map AJ can then be written in the form ( [l7| ) 
with 


a(t) = |GV )| 2 , c(t) = G(t ), (30) 


where 


7 (t) = — 25R 


C(£) 

o(ty 


(31) 


the function G(t) depending on the reservoir spectral 
density. For a qubit interacting resonantly with a leaky 
cavity, the spectral density has a Lorentzian shape, i.e, 


J(lu) — 7mA 2 /27t[(cj — cjo ) 2 + A 2 ]. 
In this case, the function G(t) takes the form 


G L (t) 


,-At/2 


cosh ( Vl- 2 ^ 


1 ( r- -A t 

. smh VI — 2 r — 

V 2 


(32) 


(33) 
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with r = 7 m/A. For the photonic band bap model (PBG) 

[35] . we have instead 

G p {t ) = 2v 1 b 1 e l3b i+ iApt + v 2 (b 2 + |5 2 |)e^ t+iApt 

- X>N[i - (34) 

3=1 

where Ap = ujq — cj e is the detuning from the band gap 
edge frequency u; e , set to equal zero as we consider only 
the resonant case, and <f> is the error function, whose 
series and asymptotic representations are given in Ref. 

[36] , and 

b 1 = {A + +A_)e i{ ^ A \ 


have plotted the incompatibility for two initial measure¬ 
ment pairs and entanglement of the initial maximally en¬ 
tangled state for the dephasing dynamics. We observe 
that the incompatibility is extremely fragile under noise. 
Indeed, entanglement decreases asymptotically towards 
zero, but incompatibility reaches zero far before the en¬ 
tanglement has vanished. Naturally, for the phase damp¬ 
ing evolution, the incompatibility of the measurements in 
the z-direction (initial pair V 2 ) is better preserved. The 
phase damping dynamics does not have a unique fixed 
point, and hence the the measurement in the z-direction 
does not evolve towards a trivial measurement. However, 
the pair V\ experiences a sudden death of incompatibility 
at an early stage of the evolution. 


62 = (A+e“ i(7r/6) — ^4_e*( 7 r / 6 ))e -i( - 7r / 4 - ) ’ 

63 = (A|_e i(7r/6) — ^ 4 _e -i ( 7r / 6 ))e*( 37r/4) , 


A+ = 


vi = 


V2 = 


Xi 


(61 - 6 2 )(6i - 63) 

b2 


fa-bjfa-h)' 

yd 3 / 2 = UJ 7 q 2 d 2 /67T€ohc 3 . 


(35) 


(37) 

(38) 

(39) 


The coefficient /3 is defined as the characteristic fre¬ 
quency, eo the Coulomb constant, d the atomic dipole 
moment, and z = Ap//3. 


VI. RESULTS AND DISCUSSION 



1 1 

4 Ap" 

1/2- 

1/3 

(36) 


2 ± 2 

1 H-— 

L 27 ^ \ 


, 

E 

0 




FIG. 1. (Color online) The evolution of entanglement en¬ 
tanglement in terms of concurrence 8 (black squares), and 
the incompatibility measure Zb=o for the dephasing evolution 
with Markovian parameters, A — 0 and Acj = 2 a. Incompat¬ 
ibility is represented by two initial pairs Vi (red circles) and 
V 2 (blue stars). 


In this section we study the dynamics of incompati¬ 
bility for the different physically relevant dynamical pro¬ 
cesses specified in the previous section. We are especially 
interested in comparing the dynamics of incompatibility 
with that of entanglement, so as to detect the fundamen¬ 
tal differences between these quantum resources. In par¬ 
ticular, we investigate the expected difference between 
the timescales of the sudden death of entanglement and 
incompatibility. 


In Fig. [2] we have plotted the incompatibility and con¬ 
currence for the Markovian amplitude damping dynam¬ 
ics. We observe the same behaviour as for the dephasing 
dynamics: a sudden death of incompatibility even with a 
fairly small amount of noise. For the amplitude damping 
case, also the initial pair V 2 exhibits a sudden death of 
incompatibility. 

B. Dynamics of Incompatibility for Non-Markovian 
Dynamics 


A. Sudden Death of Incompatibility for Markovian 
Dynamics 

Let us first study the dynamics of incompatibility for 
the case of Markovian dynamics, using the models pre¬ 
sented above. For the models considered here we are 
in the Markovian regime as long as the entanglement 
is a monotonically decreasing function. In Fig. [l] we 


As we saw in the previous section, incompatibility 
is very fragile under Markovian noise. However, quan¬ 
tum features of a system may often be recovered if the 
environment is tuned such that the dynamics exhibits 
non-Markovian behaviour [33] . Hence, in this section 
we study the dynamics of incompatibility under non- 
Markovian noise, and analyse to what extent the quan¬ 
tumness of measurements can be recovered in this case. 
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FIG. 2. (Color online) The evolution of entanglement in 
terms of concurrence £ (black squares), and the incompatibil¬ 
ity measure Tb =o for the amplitude damping dynamics with 
Lorentzian spectrum and r — 0.2. Incompatibility is repre¬ 
sented by two initial pairs Vi (red circles) and V 2 (blue stars). 

In Fig. [3] we have plotted the evolution of the incompat¬ 
ibility and entanglement for the phase damping dynamics 
in the non-Markovian case. We observe that even with 
an environment exhibiting considerable non-Markovian 
character, the incompatibility cannot be very well recov¬ 
ered (Fig. [ 3 ] i)). Naturally, for the phase damping case, 
the measurement pair with the z-direction measurement 
is again less fragile. We observe that even though a sig¬ 
nificant portion of entanglement can be recovered, the 
incompatibility remains absent. When non-Markovianity 
is increased to the extent that entanglement can be fully 
recovered, also the incompatibility appears to reach again 
its maximum value. 

For the amplitude damping evolution (Fig. [4| we ob¬ 
serve that even with highly non-Markovian dynamics, the 
incompatibility remains zero or at relatively small values, 
and cannot be fully recovered. 

We then proceed to consider a highly engineered envi¬ 
ronment, where a photonics band gap (PBG) tunes the 
dynamics such that as the non-Markovianity increases, 
and the initial entanglement can be very well protected. 
In Fig. © we plot the entanglement and incompatibility 
dynamics for this model, observing that even in the case 
where entanglement can be well protected (Fig. |5]ii)) the 
incompatibility experiences a sudden death without re¬ 
covery. Naturally, if the entanglement is fully maintained 
also the incompatibility persists. 


VII. CONCLUSIONS 

Since the discovery that quantum phenomena can be 
harnessed as resources for computation and communica¬ 
tion, a considerable effort has been made towards un¬ 
derstanding how to protect quantum systems from the 



FIG. 3. (Color online) The evolution of the entanglement in 
terms of concurrence £ (black squares), and the incompatibil¬ 
ity measure Zb=o for the phase damping dynamics with A = 1 
and Auj equal to i) 4a, ii) 8a, iii) 12cr and iv) 16cr. Incompat¬ 
ibility is represented by two initial pairs Vi (red circles) and 
V 2 (blue stars). 



A t A t 


FIG. 4. (Color online) The evolution of the entanglement in 
terms of concurrence £ (black squares), and the incompatibil¬ 
ity measure Tb=o, for the amplitude damping dynamics with 
Lorentzian spectrum and r = 50, 80, 100 and 150. Incompat¬ 
ibility is represented by two initial pairs Vi (red circles) and 
V 2 (blue stars). 

effects of noise. The basic understanding is that quan¬ 
tum properties tend to get destroyed as the system in¬ 
teracts with an environment, and hence extensive quan¬ 
titative studies on the behaviour of quantum resources 
such as entanglement and discord under noisy dynamics 
have been undertaken. However, a definite understand¬ 
ing of the essence of quantum resources remains elusive, 
as the relevant quantifications crucially depend on the 























FIG. 5. (Color online) The evolution of the entanglement in 
terms of concurrence 8 (black squares), and the incompatibil¬ 
ity measure Xb= o, for the amplitude damped dynamics with 
PBG spectral density and z — 1, 0, —5 and —10. Incompat¬ 
ibly is represented the two initial pairs Vi (red circles) and 
V 2 (blue stars). 


application at hand. In contrast to the usual viewpoint 
where quantum resources are associated exclusively to 
states, we emphasise that non-classicality may also be 
attributed to measurements used to extract information 
on the quantum system. This becomes especially rel¬ 
evant in practical scenarios where the set of available 
measurements is strongly limited by experimental con¬ 
straints. Thus, it is crucial to study how the quantum 
properties of the measurements are influenced by noise. 

In this paper, we investigated the evolution of incom¬ 
patibility of a pair of quantum measurements under two 
commonly studied dynamical maps, dephasing (in the z- 
direction) and amplitude damping. While these hardly 
exhaust the collection of physically relevant two-level dy¬ 
namics, they already represent different basic aspects of 
lossy quantum dynamics, with dephasing consisting of 
pure decoherence, and amplitude damping including also 
dissipation. More general treatment will be a topic of 
future work. 

Our starting point was the general result on incom¬ 
patibility breaking channels he which says that sudden 
death of initially maximal entanglement always implies 
sudden death of all incompatibility, while the converse is 
not true. We investigated this phenomenon in the spe¬ 
cific physically relevant settings. In particular, we found 
that, indeed, even in the absence of sudden death of en¬ 
tanglement, a sudden death of incompatibility often oc¬ 
curs. Incompatibility is known to be a crucial resource for 
many quantum protocols for which entanglement is not 
sufficient. Thus, it is possible that even in the absence 
of entanglement sudden death the relevant quantum re¬ 
source may have been already destroyed. 

We also studied how non-Markovian dynamics al¬ 


lows one to recover incompatibility. We found that in 
many cases where a significant “portion” of entangle¬ 
ment can be recovered, incompatibility cannot. Thus, 
non-Markovian dynamics may not be able to recover all 
quantum resources, even if entanglement is recovered. 
We believe that this study will be useful for developing a 
more resource-oriented view of non-Markovian dynamics. 

In contrast to concurrence and other entanglement 
measures, our choice of the incompatibility measure de¬ 
pends on the initial measurements; we demonstrated the 
dependence by using two different initial measurement 
pairs. This basis dependence is actually motivated by the 
fact that the relevant applications mentioned above (i.e. 
CHSH-Bell and steering scenarios) specifically involve a 
restricted selection of measurements. In this sense, evo¬ 
lution of incompatibility should be seen as a purpose- 
oriented view to decoherence. For instance, we could 
quantify noise-robustness of a steering-based key distri¬ 
bution protocol by incompatibility. Concerning our spe¬ 
cific choice of the incompatibility measure he there are 
alternatives, such as “incompatibility weight” [18] based 
on the steering quantification in [7j. The resulting in¬ 
compatibility measure has recently been compared with 
ours in [9] (see especially Fig. 2 therein); while there are 
strong quantitative differences, the qualitative behaviour 
is relatively similar except for special cases. 

Finally, we wish to point out that our work on the dy¬ 
namics of incompatibility in open systems has a natural 
application in quantum optimal control, namely finding 
optimal measurement resources for, say, noisy steering 
by using the incompatibility quantification as a figure of 
merit. Such a study has recently been carried out by one 
of the present authors ISZI- 
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